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I discuss NRQCD and, in particular, the NRQCD factorization formalism for quarkonium 
production and decay. I also summarize the current status of the comparison between 
the predictions of NRQCD factorization and experimental measurements. 

1. Nonrelativistic QCD (NRQCD) 

In a heavy quarkonium, a bound state of a heavy quark and antiquark, there are 
many important momentum scales. These include m, the heavy-quark mass; mv, the 
typical heavy-quark momentum in the quarkonium rest frame; and mv 2 , the typical 
heavy-quark rest-frame kinetic energy and binding energy. Here v is the typical 
heavy-quark velocity in the quarkonium rest frame. For charmonium, v 2 w 0.3, 
while for bottomonium, v 2 « 0.1. 

In theoretical analyses, it is useful to treat the physics at each of these scales 
separately. Owing to asymptotic freedom, interactions at the scale m can be treated 
perturbatively. Approximate symmetries (e.g. heavy-quark spin symmetry) can be 
exploited at some scales. Often, analytic calculations simplify when they involve 
only one scale at a time. Lattice calculations can encompass only a limited range of 
scales, and so become more tractable after scale separation. 

Effective field theories provide a convenient way to separate scales. An effective 
theory describes the low-momentum degrees of freedom in the original theory. It 
is constructed by integrating out the high-momentum degrees of freedom in that 
theory. Nonrelativistic QCD (NRQCD) is an effective field t heory that separates 
scales of order m and higher from the other scales in QCD! ' ™< NRQCD has a 
UV cutoff A ~ m. For processes with p < A, NRQCD reproduces QCD. Processes 
with p > A are not manifest in NRQCD, but they affect the coefficients of local 
interactions. A plays the role of a factorization scale between the hard and soft 
physics. 

At leading order in v, the NRQCD action is just the Schrodinger-Pauli action 
for a heavy quark and a heavy antiquark. In order to reproduce QCD completely, 
one would need an infinite number of interactions of all orders in v. However, in 
practice one works to a given precision in v. 

NRQCD and NRQED (the corresponding effective theory for QED) have had 
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a number of well-established successes. These include calculations of the properties 
of QED bound states, such as positronium and muonium, lattice calculations of 
quarkonium spectra, and lattice calculations of a s . However, here we will focus on 
recent applications to quarkonium decay and production. 

2. Quarkonium Decay and Production 

In heavy-quarkonium decays and hard-scattering production, large scales appear: 
Both the heavy-quark mass m and the quarkonium transverse momentum px are 
much larger than Aqcd- The hope is that NRQCD would allow one to separate these 
short-distance, perturbative scales from the long-distance quarkonium dynamics. 

In the case of quarkonium decays, convincing arguments have been given that 
the short-distance decay process can be represented in NRQCD by point-like four- 
fermion interactions?! This result leads to the factorization formula 



The F n are "short-distance coefficients," which are determined by perturbative 
matching of amplitudes between QCD and NRQCD. The factors (H\O n (A)\H) 
are inherently nonperturbative matrix elements of four-fermion operators in the 
quarkonium state. The sum in Eq. is actually an expansion in powers of v and 
is truncated at some finite order. A similar factorization formula has been conjec- 
tured to hold for inclusive quarkonium production at large py3 The production 
matrix elements are the crossed versions of the decay matrix elements. Only the 
color-singlet production and decay matrix elements are simply related. Nayak, Qiu, 
and Sterman have found that, in next-to-next-to-leading order in a s , the NRQCD 
production matrix elements must be modified by the inclusion of eikonal lines in 
order to allow one to factor all of the IR-divergent contributions out of the short- 
distance coefficients This finding does not affect existing phenomenology, which 
is at the tree level or one-loop level. However, it raises important issues with regard 
to the validity of the NRQCD factorization formula at all orders in a s . 

The NRQCD factorization formalism gains much of its predictive power from 
the fact that the nonperturbative matrix elements are universal, i.e., process inde- 
pendent. Although some decay matrix elements have been computed on the lattice,^ 
in general, the matrix elements must be extracted phenomenologically. The consis- 
tency of the phenomenological matrix elements from process to process is a key test 
of the NRQCD factorization formalism. 

An important feature of the NRQCD factorization formalism is that quarkonium 
decay and production occur through color-octet, as well as color-singlet, QQ states. 
If one drops all of the color-octet contributions, then the result is the color-singlet 
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Fig. 1. The J/tp cross section as a function of py. The data points are from the CDF 
measurement!^! The solid curve is the NRQCD factorization fit to the data given in Ref. 1101 
The other curves give various contributions to the NRQCD factorization fit. From Ref. 10 



3. Comparison of NRQCD factorization with experiment 
3.1. P-wave decays 

A global fit of the NRQCD predictions for the inclusive decay rates of P-wave 
quarkonia (Ref. is in good overall agreement with the data and yields values for 
the NRQCD matrix elements that are in good agreement with those from lattice 
calculations (Ref. EJ- More recently, a number of next-to- leading-order NRQCD 
predictions for decay rates of P-wave quarkonia have been verified by more precise 
experimental measurements 



3.2. Quarkonium production at the Tevatron 

The CDF data for J/tp production at the Tevatron are shown in Fig.Q along with 
the NRQCD factorization result. As can be seen, the color-singlet contributions, 
whose normalizations are reasonably well known from decay processes, are smaller 
than the data by more than an order of magnitude. The color-octet contributions 
bring the theory into good agreement with the data. However, it should be remem- 
bered that the color-octet matrix elements are obtained by fitting to the Tevatron 
data. The shape of the data is consistent with NRQCD factorization, but there is a 
good deal of freedom to change the predicted shape by adjusting the relative values 
of the color-octet matrix elements. The Tevatron data for ip' and T production are 
also fit well by the NRQCD factorization expressions. In subsections that follow, 
we discuss more stringent tests of NRQCD factorization, in which the values of the 
color-octet matrix elements that were obtained in the fits to the Tevatron data are 
used to make predictions for other processes. 
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Fig. 2. The cross section for 77 — * J /if) 4- X at LEP. The upper set of curves is the NRQCD 
factorization prediction, and the lower set of curves is t he c olor-singlet-model prediction. The data 
points are from the Delphi Collaboration measurement^^! From Ref. 1111 

3.3. 77 -> J/ip + X at LEP 

As can be seen in Fig. [21 in the case of J/ip production in 77 collisions at LEP, 
comparison of theory with Delphi data clearly favors NRQCD factorization over 
the color-singlet modelPl The theoretical uncertainties that are shown arise from 
uncertainties in the color-octet matrix elements and in the choices of renormalization 
and factorization scales. The uncertainties in the color-octet matrix elements are 
quite large because different linear combinations of matrix elements appear in J/ip 
production at LEP than in J /if) production at the Tevatron. 

3.4. Quarkonium production in DIS at HERA 

The leading-order NRQCD factorization and CSM predictions (Ref. IKS[1 for the 
J/ip inclusive production cross sections da/dp\ and da/dQ 2 in ep deep-inelastic 
scattering (DIS) are in good agreement with the HI data (Ref. I14H . Those data lie 
well above the CSM prediction. The cross section da/dz, which is differential in the 
energy fraction (inelasticity) z, is not fit well by either the NRQCD factorization 
or CSM predictions. The poor fit of the NRQCD factorization prediction may be a 
consequence of the breakdown of both the v expansion and the a s expansion near 
z = 1. This phenomenon will be discussed in conjunction with inelastic quarkonium 
photoproduction at HERA below. The ZEUS data for da/dQ 2 agree less well with 
the NRQCD prediction than the HI data, but they have larger error barsP^ 

3.5. Polarization of quarkonia at the Tevatron 

The polarizations of quarkonia produced at large pr at the Tevatron provide po- 
tentially definitive tests of the color-octet mechanism. J ftp production at large-px 
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Fig. 3. J/ip polarization at the Tevatron. a = 1 corresponds to 100% transverse polarization; 
a = — 1 corresponds to 100% longitudinal polarization. The band is the total NRQCD factorization 
prediction. The other curves give the contributions from feeddown from higher charmonium states. 
The data points are from the CDF measurement ^Hl From Ref. EE] 

is expected to be dominated by gluon fragmentation into a color-octet QQ pair. 
This mechanism leads to transversely polarized j/^'gllSIlII The NRQCD factor- 
ization prediction for the J/ip polarization as a function of pt (Ref. I18fl is shown, 
along with the CDF data, in Fig. [3J The observed J/ip polarization is generally 
smaller than the prediction and seems to trend in the wrong direction, decreasing 
with increasing px- However, the experimental error bars are large, and only the 
last data point truly disagrees with the prediction. There arc large uncertainties 
in the theoretical predictions that arise from uncertainties in the NRQCD matrix 
elements, which are indicated in the prediction band in Fig. [3J There are also large 
corrections of higher order in a, and v to the quarkonium production rate, some 
of which have been calculated!^ To first approximation, such corrections merely 
change the normalizations of the fitted color-octet matrix elements without chang- 
ing the polarization strongly. It has been suggested that nonperturbative spin-flip 
processes, which are suppressed as v and are not taken into account in present 
calculations, might be important!*^ However, a recent lattice calculation suggests 
that this is not the case!*™ 

3.6. Inelastic photoproduction at HERA 

Theoretical calculations of the cross section for inelastic photoproduction of quarko- 
nium at HERA have been carried out in the NRQCD factorization formalism by 
several groups^! The compilation of predictions from Ref. 1101 and the HI and Zeus 
data are shown in Fig. 0J plotted as function of the energy fraction z. As can be 
seen, the color-octet contribution is poorly determined, owing to large uncertainties 
in the color-octet matrix elements. Even so, there is little room for a color-octet 
contribution. Furthermore, as is shown in Fig. 0] corrections of next-to-leading or- 
der in a s (NLO) increase the color-singlet piece by about a factor of two at large 
z. The color-singlet piece is then, by itself, in good agreement with the data!^l2] 
The data differential in pr are also compatible with NLO color-singlet production 
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Fig. 4. The rate for inelastic quarkonium photo-production at HERA as a function of the energy 
fraction z. In the left-hand figure, the curves give leading-order (LO) total, direct, and resolved 
color-octet (CO) and color-singlet (CS) contributions. The band shows the uncertainty in the total 
contribution that arises from the uncertainties in the color-octet matrix elements. In the right- 
hand figure, the curves show the prediction in NLO of the color-singlet model. The band shows 
the uncertainties that arise from a s and m. The data points are from the HI and Zeus results of 
Refs . I2T1 and From Ref . ITUI 

alone at large p^P^Ilt should be noted, though, that there are large uncertainties in 
the NLO color-singlet contribution, which arise primarily from uncertainties in m c 
and a s . The true color-singlet contribution could be lower than the central value 
by about a factor of two, leaving more room for a color-octet contribution. 

Near z = 1, the leading-order color-octet contribution grows rapidly, in apparent 
disagreement with the data. However, in this region, soft-gluon emission leads to 
large logarithms of 1 — z and also to large corrections of higher order in v, both of 
which must be resummed. The resummation of the corrections of higher order in 
v leads to a nonperturbative "shape function. Both the shape function and the 
resummed logarithmic corrections significantly smear out the color-octet contribu- 
tion near z = 1 and may lead to a considerable improvement in the agreement of 
the NRQCD factorization predictions with the data. 

3.7. Double cc production at Belle 

For the exclusive double charmonium process e + e~ — -> J/ip + i] c , the cross section 
times the branching ratio into at least two charged tracks has been measured by the 
Belle collaboration to be 25.6±2.8±3.4 fb (Ref. 26) and by the BaBar collaboration 
to be 17.6 ± 2.8^21 fb (Ref. EH- In contrast, leadin g-orde r NRQCD factorization 
calculations predict a cross section of 3.78 ± 1.26 A similar disagreement 

between NRQCD factorization and experiment holds for production of \c0 and 
r] c (2S) mesons in conjunction with a J/ip meson. A recent calculation of corrections 
of next-to-leading order in a s leads to an enhancement of the theoretical prediction 
of about a factor 1.8 (Ref. I30[) . Calculations in the light-cone formalism are in 
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reasonably good agreement with the dataP^^l The light-cone formalism takes into 
account effects from the relative motion of the Q and Q in the quarkonium, which 
is neglected at leading order in NRQCD. However, it is not clear that the model 
wave functions that are used in these calculations are good approximations to the 
true quarkonium light-cone wave functions. 

There are also Belle results on inclusive doublc-charmonium production. For 
the ratio Rj/^p = a(e + e~ — ► J/ip + cc)/<r(e + e~ — > J/tp + X), the most recent Belle 
analysis yields R J/4 , = 0.82±0.15±0.14, with R Jf ^ > 0.48 (90% confidence level) 
Predictions based on NRQCD factorization give Rj/^p ~ 0.1 fRef. I17jl . In the case 
of the absolute cross section for J/tp + cc production, the Belle result of 0.6-1.1 pb 
(Ref. 154"}) disagrees with the prediction of 0.10-0.15 pb (Ref. ITY|l by almost an 
order of magnitude. This prediction is based only on the color-singlet contribution. 
However, corrections of higher order in v, including color-octet contributions, are 
not expected to be large. Neither are corrections of higher order in a s . It is difficult 
to see how any perturbative calculation could give a value for Rj/^ as large as 80%. 

4. Summary 

The effective field theory NRQCD is a convenient formalism for separating physics 
at the scale of the heavy-quark mass from physics at the scale of quarkonium bound- 
state dynamics. The NRQCD factorization approach provides a systematic method 
for calculating quarkonium decay and production rates as a double expansion in 
powers of a s and v. NRQCD factorization for production rates relies upon hard- 
scattering factorization and has not yet been established. The NRQCD factorization 
approach has enjoyed a number of successes in inclusive P-wave quarkonium decays, 
quarkonium production at the Tevatron, 77 — > J/ip + X at LEP, and quarkonium 
production in DIS at HERA. Other processes, including production of polarized 
quarkonium at the Tevatron, inelastic quarkonium photoproduction at HERA, and 
double cc production at Belle and BaBar are more problematic and point to the fact 
that our theoretical understanding of quarkonium production is still incomplete. 
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